F/G  15/5 


AD-A080  952  ARMY  INVENTORY  RESEARCH  OFFICE  PHILADELPHIA  PA 
A  NOTE  ON  INITIAL  FILL  RATE* (U) 

FEB  BO  A  J  KAPLAN 

UNCLASSIFIFD  IR0-TR-80-1  NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BURLAU  OF  STANDAI  OS  1963  A 


htemin  Hi  data  tarn 
■■y  W  subnet  it  claaga,  tkla  4aei 
ualtas  «•  ttaw4. 


SECURITY  CLASSIFICATION  OF  THIS  PAOE  (Whan  Data  Enlarad) 


REPORT  DOCUMENTATION  PAGE 


I.  REPORT  NUMBER 


%Arc  READ  INSTRUCTIONS 

_ BEFORE  COMPLETING  FORM 

2.  GOVT  ACCESSION  NO.  »■  RECIPIENT'S  CATALOG  HUMaER 


I*.  TITLE  (and  Sublllla) 


S.  TYPE  OF  REPORT  A  PERIOD  COVEREO 


Technical 


A  NOTE  ON  ^INITIAL^ILL  RATE 


F*  AUTHORfaJ 


I  A.  CONTRACT  OR  GRANT  NUMSENfaJ 


J}Ji N  JJ  KAPLAN  1 _ 

pZVLUfflnim  l.llll  llu/lTION  NAME  AND  AOORESS  ,  10.  PROGR AM^EMEI^T.JMMy ECT,  TASK 

US  Army  Inventory  Research  Office/ ALMC 
Room  800,  US  Custom  House,  2nd  &  Chestnut  Sts. 

Philadelphia,  PA  19106  _  _ _ 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  "l  ii  IILI  IJR  I  U  ATI  1 

US  Army  Materiel  Development  &  Readiness  Command  /  //  FebMBapaM8y  / 

5001  Eisenhower  Avenue  '''--♦4:  *(JtoBfRol-l'pAGEl 

Alexandria,  VA  22333 _ _ 14 _ 

U  MONITORING  AGENCY  NAME  A  ADDRESSfi/  dlllarant  from  Controlling  Ottlca)  IS.  SECURITY  CLASS,  (of  thla  report) 


/Mr 


UNCLASSIFIED _ 

15*.  declassification/downgrading 
SCHEDULE 


1 16.  DISTRIBUTION  STATEMENT  (of  fhfl 


Approved  for  Public  Release;  Distribution  Unlimited 


1 17.  DISTRIBUTION  STATEMENT  (of  tho  mbatrocl  antarad  In  Block  20,  li  dltloront  from  Roport) 


/±1 


is.  supplementary  notes  /  _  _  — — — - - - / 

Information  and  dats^contalned  in  this  document  are  based  on  input  avail¬ 
able  at  the  time  of  preparation.  Because  the  results  may  be  subject  to  change, 
this  document  should  not  be  construed  to  represent  the  official  position  of 
the  US  Army  Materiel  Development  &  Readiness  Command  unless  so  stated. _ 

IS.  KEY  WORDS  (Conllnua  an  tararaa  alda  II  racaaaarr  and  Idanlllr  by  block  numb at) 

E0Q 

Reorder  Points 
Performance 


20.  ABSTRACT  (Caotbrua  am  rararaa  NA  H  naaaaaaay  and  Identity  by  block  mmbar) 

*^A  common  measure  of  customer  service  is  initial  fill  rate,  the  percent  of 
demand  satisfied  without  delay  from  stock  on  hand.  An  exact  expression  for 
initial  fill  rate  when  demand  is  Compound  Poisson  has  previously  been  developed. 
This  note  presents  an  approximation  which  is  simple  to  program  and  relatively 
quick  to  run.  Its  accuracy  is  distribution  dependent.  For  Negative  Bimmlal 
distribution  of  demand,  its  accuracy  is  evaluated  and  found  to  be  good.  A 
common  alternative  approach  to  calculate  Initial  fill  is  shown  to  work  very 
badly  for  one  class  of  items  _ _ _ 

DO  /jZTn  U73  TwrwAQMiioyiiii^woLETE  UNCLASSIFIED  ^7  *)  ^ 

SECURITY  CLASSIFICATION  OF  TRIE  PAGE  (man  bataTntarad) 


TABLE  OF  CONTENTS 


TABLE  OP  CONTENTS .  1 

1.  SUMMARY .  2 

2.  APPROXIMATION .  2 

3.  COMPUTATIONAL  CONSIDERATIONS .  4 

4.  EVALUATION .  5 

5.  CONCLUDING  REMARKS . 8 

REFERENCES..... .  9 

APPENDIX . 10 

DISTRIBUTION . . .  12 


I 


DTIC 

s 

I  E 


ELECTE 

FEB  2  2  1980 1 

B 


I  ACCESSION  for  L 

NTIS 

DOC 

UNANMOUN 

WSTtnCM 

Whitt  Sect  ion 

Buff  StetioR  □ 

ecu  □ 

inn 

91 

OBDtfBtU 

Dnt.  A 

IW/AVMLMUTY  HO 

/AIL.  and/or  SPECIAL 

jt\ 

I 


A  NOTE  ON  INITIAL  FILL  RATE 


1.  Summary 

A  common  measure  of  customer  service  is  initial  fill  rate,  the  percent 
of  demand  satisfied  without  delay  from  stock  on  hand.  Silver  has  developed 
an  exact  expression  for  initial  fill  rate  when  demand  is  Compound  Poisson. 
This  note  presents  an  approximation  which  ia  far  simpler  to  program,  and 
much  less  time  consuming  to  run.  Its  accuracy  is  distribution  dependent. 

For  Negative  Binomial  distribution  of  demand,  its  accuracy  is  evaluated 
and  found  to  be  good.  The  Negative  Binomial  was  chosen  because  of  its  use 
in  a  standard  US  Army  inventory  system  (3) ,  and  its  general  suitability 
for  modelling  demand. 

Other  assumptions  are  that  unfilled  demand  is  backlogged,  that  the 
Negative  Binomial  parameters  are  known,  and  that  a  continuous  review  reorder 
point,  order-up-to-level  policy  is  used.  Because  transaction  size  may  exceed 
one,  inventory  at  the  time  of  an  order  may  be  below  the  reorder  point.  The 
difference  is  termed  the  deficit. 

Notation: 

R  -  reorder  point 
M  -  order-up-to-level 

Q  -  reorder  quantity  for  zero  deficit  (M-R) 

Z  -  deficit 
X  -  lead  time  demand 
S  -  transaction  size 

f(x),  t(z),  (h(s)  are  probability  mass  functions  for  X,  Z,  S. 

IF(r,q)  =  initial  fill  rate  assuming  every  deficit  is  zero,  the 
reorder  point  is  r,  and  the  reorder  quantity  is  q. 

G(r,q)  =  IF(r,q)  •  q 
Z1  =  Ex  (zlz  <  R) 

2.  Approximation 

Full  understanding  of  this  section  requires  a  review  of  Silver's  (6) 
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paper.  From  Silver's  paper  (equation  top  of  p.96),  it  follows 


(rfq-x)f(x) 


By  using  the  function  G(r,q),  Silver's  equation  (11)  can  be  reveallngly 
rewritten  as: 


E  t(z)G(R-z,Q*-*)  +  pr(z  >  R)G(0,R+Q) 
z»o 


Initial  Fill  Rate 


Q+E(Z) 


t(z)G(R-z,Q+z)  +  pr(z  >  R)G(O.R-K)) 


Q+E(Z) 


The  first  step  in  the  approximation  is  to  set 


t(z)G(R-z,Q+z)  ~Pr(z  <  R)G(R-Z  ) 


The  second  step  in  the  approximation  is  to  represent  the  distribution 
of  the  number  of  units  demanded  per  transaction  by  the  geometric  distribution 
For  the  geometric  it  is  known  (cf  Silver,  Ho,  Deemer)  that  the  distribution 
of  the  deficit  is  the  same  as  the  distribution  of  transaction  size  plus  1; 
l.e.  if  W  is  the  amount  by  which  inventory  exceeds  R, 


[h(s+w-l)|s  >  v-1 ]  ■  h(s)  ■  t(s-l) 
Defining  the  geometric  function: 

/  _  \  8“1  _ j i 


I  first  saw  these  results  in  notes  by  Richard  Urbach 
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One  problem  which  does  arise  Is  that  R-Z^ ,  Q+Z^  need  not  be  Integer 
while  the  negative  Binomial  Is  defined  only  for  integer  values.  Our  solu¬ 
tion  was  to  interpolate:  e.g.  if  "fr"  is  a  fractional  value 

(8)  f(*ffr;  w,u.v)  -  (fr)  f(xfl;w,u,v)  +  (1-fr)  f(x;w,u,v) 

Camp  Paulson  Transformation.  The  Negative  Binomial  Cumulative 

Distribution  Function  may  be  approximated  by  the  Normal,  provided  the 

value  at  which  the  c.d.f.  is  evaluated  is  transformed  in  the  manner  devised 
* 

by  Camp  Paulson.  (In  other  words  if  we  wish  to  know  the  probability  demand 

i 

is  <  X,  we  transform  X  to  X  and  find  the  probability  a  standard  normal 

~  »  (2) 
is  <  X  ) .  Furthermore  , 

B  B-l 

(9)  l  x  f (x;w,u,v)  «  ^  i  f(x;w+l,u,v) 

x-o  X“0 

Letting  B  -  M,  and  then  R,  we  see  that  G(r,q)  can  be  expressed  as  a  function 
of  Negative  Binomial  c.d.f. 's,  and  then  evaluated  using  the  Camp  Paulson 
approximation.  The  Camp  Paulson  approach  is  useful  when  lead  time  demand 
can  be  high,  and  incidentally  eliminates  any  need  to  interpolate  to  treat 
non-integer  values  as  in  equation  (8) . 

Computational  Comparison  of  Approaches.  Even  if  the  deficit  is  ignored 
G(r,q)  must  be  computed.  The  approximate  algorithm  requires  evaluation  of 
G(r,q)  at  two  different  points  and  not  much  extra  programming.  As  compared 
to  the  exact  approach,  it  eliminates  the  need  to  compute  the  deficit  dis¬ 
tribution,  and  to  evaluate  G(r,q)  at  a  possibly  large  number  of  different 
points. 

4.  Evaluation 

The  evaluation  compared: 

"Exact"  Silver's  formula 

"Approx"  Approximation  as  given  in  equation  (5) 

See  reference  (1)  or  (2) .  In  (1) ,  the  minus  sign  before  2nd  term  in 
equation  (7)  definition  of  z,  should  read  plus. 
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"Approx/CP"  Approximation  using  Camp  Paulson 

"Ignore"  Use  IF(R,Q),  ignoring  the  impact  of  the  deficit. 

"Ignore"  might  also  have  been  characterized  as  "Poisson-Extended",  as 
an  alternative  derivation  of  "Ignore"  is: 

a.  Write  down  the  expression  for  initial  fill  under  Poisson  demand. 

b.  Replace  the  Poisson  c.d.f.  with  the  Compound  Poisson  c.d.f.  (or 
perhaps  some  other  such  as  Normal  or  Laplace) . 

The  equivalence  of  the  result  from  step  (b)  and  equation  (1)  can  be 
shown  with  some  algebra  and  calculus  (see  Appendix) .  As  "Poisson  Extended" , 
"Ignore"  is  commonly  incorporated  into  inventory  models;  e.g.,  see  (2),  (5). 

Inspection  of  the  tabulated  results  of  the  evaluation  confirms  that 
when  S  is  large  relative  to  lead  time  demand,  the  deficit  distribution 
must  be  considered.  Silver  calls  such  items  "erratic  demand"  items. 

"Approx"  and  "Approx/CP"  work  equally  well.  However,  their  error  does 
not  decrease  as  the  deficit  distribution  becomes  less  important  -  focus 
on  the  case  where  X  -  10,  S  ■  2  -  even  though  in  the  limit  as  S  -►  1,  "Approx" 
and  "Ignore"  must  converge  to  "Exact". 

The  impact  of  equation  (3),  taken  by  itself,  is  to  produce  results 
between  "Exact"  and  "Ignore".  When  more  of  the  deficit  distribution  is 
larger  than  R,  equation  (3)  has  less  Impact,  so  with  the  exception  of  the 
(10,2)  cases,  approx  does  better  for  R  ■  X  than  for  R  =  2X. 

The  other  aspect  of  the  approximation  is  use  of  the  Geometric  to 
characterize  transaction  sizes  Instead  of  the  Logarithmic  distribution 
implied  by  use  of  the  Negative  Binomial  (4)  and  incorporated  into  "Exact". 
Additional  runs  on  the  (10,2)  cases  showed  this  was  not  a  problem; 
in  fact,  the  impact  of  the  Geometric  by  Itself  was  to  significantly  under¬ 
state  initial  fill. 

Still  other  runs  not  shown  verified  that  Approx  worked  well  for 
expected  demand  and  order  size  of  (25,2),  (25,5),  (25,25),  (50,2),  (50,5), 
(50,25). 

*Even  for  the  (10,2)  case,  if  error  is  measured  as  a  !  of  unfilled  demand 
(e.g.  . 9%/<1002— 94 . 1%) ,  Approx  does  better  when  R  ■  X. 
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COMPARISON  OF  APPROACHES 


X 

S 

R 

Q 

"Exact" 

"Approx" 

"Approx/CP" 

"Ignore" 

2 

2 

X 

2S 

70. 2X 

+0.6% 

-0.6% 

+11.8% 

X 

5S 

84. 2X 

+0.4% 

-0.3% 

+  6.6% 

2X 

2S 

83.9% 

+1.5% 

+1.8% 

+  7.4% 

2X 

5S 

91.5% 

+0.8% 

+0.9% 

+  4.1% 

2 

4 

X 

2S 

64.1% 

-0.4% 

+0.0% 

+24.2% 

X 

5S 

79.8% 

+0.3% 

+0.3% 

+14.0% 

2X 

2S 

72.5% 

+0.6% 

+1.2% 

+19.4% 

2X 

5S 

84.5% 

+0.7  % 

+0.8% 

+11.1% 

10 

2 

X 

2S 

62.2% 

+2.1% 

+2.3% 

+  6.1% 

X 

5S 

75.5% 

+1.0% 

+1.1% 

+  4.1% 

2X 

2S 

94.1% 

+0.9% 

+0.8% 

+  1.7% 

2X 

5S 

96.5% 

+0.4% 

+0.4% 

+  1.0% 

10 

5 

X 

2S 

60.1% 

+1.2% 

+1.9% 

+16.3% 

X 

5S 

74.4% 

+1.1% 

+1.3% 

+10.9% 

2X 

2S 

79.6% 

+3.0% 

+3.1% 

+  9.9% 

2X 

5S 

87.0% 

+1.9% 

+1.8% 

+  6.4% 

10 

10 

X 

2S 

56.8% 

-0.7% 

+0.1% 

+26.7% 

X 

5S 

72.5% 

+0.4% 

+0.7% 

+17.4% 

2X 

2S 

69.4% 

+1.0% 

+0.6% 

+20.2% 

2X 

55 

80.3% 

+1.4% 

+1.6% 

+13.2% 
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5. 


Concluding  Remarks 


The  approximation  described  in  this  note  is  not  fully  satisfying  because 
it  depends  in  part  on  the  ability  of  the  geometric  to  characterize  trans¬ 
action  size.  At  the  same  time  there  are  many  items  for  which  it  is  useful, 
and  for  which  ignoring  the  deficit  distribution  is  unthinkable. 

In  80 me  applications  the  variability  of  lead  time  demand  is  Increased 
as  an  approximate  way  of  coping  with  lead  time  variability  or  uncertainty 
about  the  parameters  of  the  lead  time  demand  distribution.  This  is  readily 
accommodated  by  the  approximation,  since  it  need  only  change  the  calculation 
of  G(r,q) . 
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APPENDIX 


Backup  to  Equation  4b 


P 


r  i 

la.  a 


l-l 


x  d  1  .  a  b+1] 

q  ‘Pdq  ’‘I 

_<b+l)_.qbJ(l-q)  +  qfl  -  q^1  . 

(1-q)2 


a-1  ,.N  b  ,  a  b+1 

a  q  -  (b+l)q  +  q  -q  _ 

P 

qa_1  [a  +  q/p]  -  qb[ (b+1)  -  q/p] 


Backup  to  Equation  (5)  , 


R 

I  s  h(s) 

a-1 _ 

R 

l  h(s) 
8-1 


(I-Hl/P)  -  qR[R+l  +  q/p]  _ 

l-qR 

(i)  -  qy  +  i3 


L 

» 
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Backup  to  Equivalence  of  "Ignore"  and  "Poisson  Extended" 


For  "Poisson  Extended" 


±  tQF(R)  +  l 


i  [QF(R)  +  I 


±  [QF  (R)  +  I 


F(R)  +  ±  I 


F(R)  +  ±  l 


R+Q 

I 

x-R+1 

F(x-l)  - 

R+Q 

X-l 

E 

E  f  (y)  ] 

x-R+2 

y-R+1 

R+Q-l 

X 

E 

e  f(y)J 

x-R+1 

y-R+1 

sunmation 

R+Q-l 

R+Q-l 

E 

E  f (y)] 

y-R+1 

x-y 

R+Q-l 

E 

y-R+1 

(M-l-y+1)  f (y) 

R+Q-l 

E 

y-R+1 

(M-y)f(y)  - 

’Ignore" 
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